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ABSTRACT 



The existence of ionized X-ray absorbing layers of gas along the line of sight to 
the nuclei of Seyfert galaxies is a well established observational fact. This material is 
systematically outflowing and shows a large range in parameters. However, its actual 
nature and dynamics are still not clear. In order to gain insights into these important 
issues we performed a literature search for papers reporting the parameters of the soft 
X-ray warm absorbers (WAs) in 35 type 1 Seyferts and compared their properties to 
those of the ultra-fast outflows (UFOs) detected in the same sample. The fraction of 
sources with WAs is >60%, consistent with previous studies. The fraction of sources 
with UFOs is >34%, >67% of which also show WAs. The large dynamic range obtained 
when considering all the absorbers together, spanning several orders of magnitude in 
ionization, column, velocity, and distance allows us, for the first time, to investigate 
general relations among them. In particular, we find significant correlations indicating 
that the closer the absorber is to the central black hole, the higher the ionization, 
column, outflow velocity and consequently the mechanical power. In all the cases, the 
absorbers continuously populate the whole parameter space, with the WAs and the 
UFOs lying always at the two ends of the distribution. These evidences strongly suggest 
that these absorbers, often considered of different types, could actually represent parts 
of a single large-scale stratified outflow observed at different locations from the black 
hole. The UFOs are likely launched from the inner accretion disc and the WAs at larger 
distances, such as the outer disc and/or torus. We argue that the observed parameters 
and correlations are, to date, consistent with both radiation pressure through Compton 
scattering and MHD processes contributing to the outflow acceleration, the latter 
playing a major role. Most of the absorbers, especially the UFOs, show a sufficiently 
high mechanical power (at least ~0.5% of the bolometric luminosity) to provide a 
significant contribution to AGN feedback and thus to the evolution of the host galaxy. 
In this regard, we find possible evidences for the interaction of the AGN wind with 
the surrounding environment on large-scales. 

Key words: accretion, accretion discs - black hole physics - galaxies: active - galax- 
ies: Seyfert - X-rays: galaxies 



1 INTRODUCTION 

The presence of ionized material along the line of sight to 
the nuclei of Seyfert galaxies has been known for a long time 
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(e.g. Halpern 1984). The material observable in absorption 
in the soft X-rays has been referred to as a warm absorber 
(WA). The limited energy resolution of previous X-ray in- 
struments allowed essentially only the detection of broad ab- 
sorption edges. In fact, observations with ASCA (Advanced 
Satellite for Cosmology and Astrophysics) detected O VII 



2 F. Tombesi et al. 



and O VIII absorption edges and established that WAs arc a 
common feature of AGNs, being present in ~50% of Seyfert 
galaxies (Reynolds 1997; George et al. 1998). Then, the ad- 
vent of the unprecedented spectral resolution of the gratings 
on board Chandra and XMM-Newton allowed, for the first 
time, the detection of discrete soft X-ray resonance absorp- 
tion and emission lines. The resulting picture of the WA 
is that of an outflow exhibiting multiple narrow absorption 
lines corresponding to different ionization states (Kaastra 
et al. 2000; Kaspi et al. 2000; Blustin et al. 2005; McKcr- 
nan et al. 2007). The values of the ionization parameter are 
typically in the range log^~0-2 erg cm and the column 
densities arc between iVH~10^" 10^^ cur^. The absorption 
lines are systematically blue-shifted, indicating outflow ve- 
locities of the WAs in the range Wout~100-1000 km/s. There 
are still significant uncertainties on the exact location of 
this material, which ranges from ~pc up to ~kpc scales, 
and it has been suggested that it might originate outside 
of the inner disc, probably at locations comparable with 
the obscuring torus (e.g. Krolik & Kriss 2001; Blustin et 
al. 2005). Depending on the actual filling and covering fac- 
tors, the mass outflow rate from the WAs can be significant 
but, given their relatively low velocities, their kinetic power 
is rather low when compared to the bolometric luminosity 
(e.g. Blustin et al. 2005; McKernan et al 2007). However, a 
recent detailed study by Crenshaw &: Kraemer (2012) found 
that the integrated power of the WAs in some Seyferts can 
actually reach the level of ~0. 1-0.5% of the bolometric lu- 
minosity, the minimum required by numerical simulations 
for AGN feedback to exert a significant impact on the host 
galaxy (e.g. Hopkins & Elvis 2010; Gaspari et al. 2011a, b, 
2012b). 

Besides WAs, highly blucshifted Fe K-shell absorption 
lines at E>7 keV have been detected in more recent years 
in the X-ray spectra of several AGNs (Chartas et al. 2002, 
2003; Pounds et al. 2003; Dadina et al. 2005; Markowitz et 
al. 2006; Braito et al. 2007; Cappi et al. 2009; Reeves et 
al. 2009; Chartas et al. 2009; Giustini et al. 2011; Gofford 
et al. 2011; Lobban et al. 2011; Dauser et al. 2012). In par- 
ticular, a uniform and systematic search for blueshifted Fe 
K absorption lines in a sample of 42 local (z ^ 0.1) Seyferts 
observed with XMM-Newton was performed by Tombesi et 
al. (2010a, hereafter Paper I). This allowed the authors to 
assess their global significance and derive a detection firac- 
tion of >40%. In order to mark an initial and somewhat 
arbitrary distinction with the classical WAs, in paper I we 
defined ultra-fast outflows (UFOs) as those highly ionized 
Fe K absorbers with a blueshifted velocity ^10,000 km/s. 
Subsequently, Tombesi et al. (2011a, hereafter Paper II) 
performed a photo-ionization modelling of these UFOs and 
derived the distribution of their main physical parame- 
ters. The outflow velocity is mildly-relativistic, in the range 
~0.03-0.3c (~10,000-100,000 km s"^), with mean value of 
~0.1c. The ionization is very high, in the range log^~3- 
6 erg s~^ cm, with a mean value of ~4.2 erg s^^ cm. The 
column densities are also large, in the interval Nh^W^^— 
10^* cm~^. These flndings are important because they sug- 
gest the presence of massive and highly ionized absorbers 
outflowing with mildly-relativistic velocities from the nuclei 
of these Seyfert galaxies. 

In a following paper, Tombesi et al. (2012a, hereafter 
Paper III) quantifled that they are observable at locations 



of sub-parsec scales from the central black hole. Their mass 
outflow rate was constrained between ~0.01-1 Mq yr^^ 
and their kinetic power was found in the range logi^K— 42- 
45 erg s"^ Thus, the UFOs are possibly directly identifiable, 
at least qualitatively, with accretion disc winds/outflows 
(King & Pounds 2003; Proga & Kallman 2004; Ohsuga et 
al. 2009; Sim et al. 2010; Pukumura et al. 2010) or the base 
of a possible weak jet (e.g. Ghisellini et al. 2004). In particu- 
lar, the kinetic power of these UFOs is systematically higher 
than the minimum required by sinnilatioiis of feedback in- 
duced by winds/outflows (e.g. Hopkins & Elvis 2010; Gas- 
pari et al. 2011a, b, 2012b). Therefore, in the long term, they 
could be able to significantly influence the bulge evolution, 
star formation, super-massive black hole growth and con- 
tribute to the establishment of the observed black hole-host 
galaxy relations, such as the Mbh-o (Ferrarese & Merritt 
2000; King 2010a; Ostriker et al. 2010; Gaspari et al. 2011a, 
b, 2012a, b). 

It is important to note that Tombesi et al. (2010b, 
2011b, 2012b) detected the presence of UFOs also in a 
small sample of (3 out of 5) radio-loud AGNs observed with 
Suzaku. Finally, we also note that similar results regarding 
the UFOs have been obtained independently by Gofford et 
al. (2012) who performed a uniform and systematic broad- 
band spectral analysis of a large sample of AGNs observed 
with Suzaku, conflrming their overall incidence and charac- 
teristics. Moreover, evidences for similar outflows are emerg- 
ing also in stellar-mass black holes (e.g.. King et al. 2012) 

Despite these significant observational developments, 
the origin and acceleration mechanism(s) of the ionized ab- 
sorbers in AGNs are still debated. In particular, radiation 
and MHD wind models have been developed and scenarios 
in which they are considered as intrinsically distinct or as 
different manifestations of the same phenomenon have been 
suggested (e.g. Konigl & Kartje 1994; Krolik & Kriss 2001; 
Elvis 2000; Blustin et al. 2005; Krongold et al. 2007; Fuku- 
muraet al. 2010; Kazanas et al. 2012; Reynolds 2012). Given 
the relevance of these outflows for the physics and energet- 
ics of AGNs and their potential significant contribution to 
feedback, it is imperative to investigate them in detail. 

In order to test these hypothesis, in this paper we will 
perform a detailed observational comparison of the WAs and 
UFOs, checking for correlations and discussing the possible 
unification of these absorbers in a single, photo-ionized and 
stratified outfiow. We focus on the sample of Seyfert galax- 
ies described in Paper I and we will use WA parameters 
collected from the literature and those of the UFOs derived 
in Paper II and Paper III, respectively. 



2 WARM ABSORBERS SELECTION 

We performed a literature search for papers reporting the 
analysis of the soft X-ray WAs in the 35 type 1 Seyferts of 
the sample discussed in Paper I. This was defined selecting 
all the Seyfert galaxies from the RXTE All-Sky Slew Survey 
Catalog (Revnivtsev et al. 2004) and cross-correlating tliem 
with the XMM-Newton Accepted Targets Catalog (as of Oc- 
tober 2008). After applying the standard filtering processes, 
we obtained a total of 42 objects with 101 good XMM- 
Newton observations and, more specifically, 35 classified as 
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type 1 and 7 as type 2 with 87 and 14 observations, respec- 
tively. 

In the literature search, we selected only the WA re- 
sults derived from the high energy resolution gratings on 
board Chandra and XMM-Newton because they allow us to 
constrain the outflow velocity, which is crucial to derive the 
mass outflow rate and the kinetic power. We limited our 
search only to the type 1 sources because, in accordance 
with the uniflcation model, the spectra of the 7 type 2s are 
affected by signiflcant neutral absorption (A'h^IO^^ cm~^) 
which hampers the detection of the WAs in the soft X-rays. 
For the UFOs, we consider the parameters reported in Paper 
II and III. 

In the following, we use the definition of the ioniza- 
tion parameter ^ = Lio-^jnr^ (Tarter et al. 1969) in which 
Z/ion is the ionizing luminosity between 1 Ryd and 1000 Ryd 
(1 Ryd = 13.6 eV), n is the number density of the absorb- 
ing material and r is the distance from the central source. 
Often, different WA components with diverse ionization, ve- 
locities and column densities are detected for each source 
and there might be time variability, especially between ob- 
servations spaced by several years. Moreover, there might 
be some intrinsic inhomogeneities in the density and ion- 
ization structure of the absorbing material. However, here 
we do not consider any subtle variations in the warm ab- 
sorber properties because we focus on deriving the global 
properties of the outflows and we refer the reader to the 
papers reported in the notes of Table 1 for more informa- 
tion. Therefore, for each source we report the various ionized 
absorption components in Table 1. If more than one paper 
reported the analysis of the same source, we averaged the 
values of the components with equivalent parameters. This 
allows us to minimize the scatter due to different analy- 
sis methods employed by different authors and the effects 
of time variability as well. The WA parameters for all the 
sources, along with their central black hole masses and aver- 
age (absorption corrected) ionizing luminosities Lion derived 
from the XMM-Newton observations reported in Paper I are 
reported in Table 1. In the subsequent correlation analysis 
we will include the points from these multiple zones of warm 
absorbers separately. 

For comparison, we use the parameters of the highly 
ionized Fe K-shell absorbers reported in Paper II. These 
were initially simply distinguished as UFOs or non-UFOs if 
their velocity was higher or lower than 10,000 km/s, respec- 
tively. However, as noted in Paper II and HI, their parame- 
ter distributions are not bimodal, but they actually show a 
roughly continuous distribution in ionization, column den- 
sity and velocity, with the UFOs lying at the more extreme 
side. This point is also confirmed by the aiiatysis of Gofford 
et al. (2012). Bearing this in mind, in the following we will 
continue to refer to them as UFOs and non-UFOs (to indi- 
cate those detected in the Fe K band with an intermediate 
ionization/velocity), but they will be correctly considered 
together in the subsequent correlation analysis in §4. 

The number of sources having papers reporting the de- 
tection of WAs is 21/35, therefore the fraction of objects 
with WAs is at least >60%, consistent with previous studies 
(Reynolds 1997; George et al. 1998; Blustin et al. 2005; McK- 
ernan et al. 2007; Winter 2010). If we consider the sources 
showing absorbers in the form of WAs, UFOs or non-UFOs 
this increases to 26/35, >74%. This suggests that the abso- 



lute majority of bright Seyfert 1 galaxies do show some form 
of ionized X-ray absorption if examined with sufficiently 
high S/N observations (Winter 2010; Paper I). The fraction 
of sources showing UFOs is 12/35 (>34%) and 8/12 (>67%) 
of these show also WAs. If we consider the UFOs and non- 
UFOs together we obtain a fraction of 16/35 (>46%) and 
consequently 11/16 (>69%) of these sources show also WAs. 
These fractions might possibly depend also on the inclina- 
tion of the flow with respect to the line of sight. Considering 
the fact that some absorbers may have not been detected 
due to low S/N, variability or simply because some sources 
have no grating observations or the WAs were not studied in 
detail, we emphasize that these fractions do represent only 
lower limits. 



3 WARM ABSORBER PARAMETERS 

We estimate the lower and upper limits of the distance, mass 
outflow rate and kinetic power of the WAs following the 
method outlined in Paper HI. An upper limit on the line of 
sight projected location can be derived from the definition 
of the ionization parameter reported in §2 and the require- 
ment that the thickness of the absorber does not exceed its 
distance to the black hole, A'h — nAr < nr (e.g. Crenshaw 
& Kraemer 2012): 

rmax = Lion/^Nu, (1) 

the material can not be farther away than this given the 
observed ionization and column. Instead, an estimate of the 
minimum distance can be derived from the radius at which 
the observed velocity corresponds to the escape velocity: 

rn,in = 2GMBH/vLf (2) 

For the calculation of the mass outflow rate we use the ex- 
pression derived by Krongold et al. (2007), which is more 
appropriate for a biconical wind-like geometry: 

Mout = f{5, cj>)7riJ.mpNHVoutr, (3) 

where f{S, 0) is a function that depends on the angle be- 
tween the line of sight to the central source and the accretion 
disc plane, S, and the angle formed by the wind with the ac- 
cretion disc, (p (see Fig. 12 of Krongold et al. 2007). Instead, 
/i = nn/rie ^ 1/1 A for Solar abundances. For a roughly 
vertical disc wind {<P~it/2) and an average line of sight an- 
gle of (5~30° for the Seyfert Is considered here (Wu & Han 
2001) we have /((5, (/>)~1.5. Full details on the derivation of 
this formula can be found in the Appendix 2 of Krongold et 
al. (2007). 

This expression for the mass outflow rate has also the 
important advantage of not relying on the estimate of the 
covering and filling factors. This is due to the fact that it 
takes into account only the net observed thickness of the 
gas, allowing for clumping in the flow. Thus, there is not the 
need to include a linear (or volume) flUing factor, since we 
are interested in estimating the net flow of mass, starting 
from the observed column density and velocity. Moreover, 
the covering factor is implicitly taken into account by the 
function f{5, 0) when calculating the area filled by the gas, 
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Table 1. Parameters of the soft X-ray WAs for the type 1 Seyferts in the sample. 
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Notes. C and X stand for grating observations with Chandra or XMM-Newton, respectively. The * and + mark the sources with detected Fe K absorbers identified 

with UFOs and non-UFOs in Paper II, respectively. ^ Kraemer et al. (2005); ^ McKernan et al. (2007); ^ Steenbrugge et al. (2005); * Matt et al. (2006); ^ Yaqoob et 

al. (2003); ® Vaughan ct al. (2004); Cardaci ct al. (2011); * Blustin ot al. (2003); ^ Svoboda ct al. (2010); Costantini et al. (2007); Markowitz et al. (2009); 

1^ Starling ct al. (2005); Gallo ot al. (2011); Lobban ct al. (2011); ^'^ Pounds & Vaughan (2011); Longinotti ct al. (2010); Laha ct al. (2011); 

Emmanoulopoulos ct al. (2011) .Timcncz-Bailon ct al. (2008); ^" Pounds ct al. (2004); Gondoin ct al. (2002); Porquct ct al. (2004); Zhang et al. (2011); 

Reeves ct al. (2001); Longinotti ct al. (2007); Ashton ct al. (2006); Hozanska ct al. (2004); Pounds ct al. (2003); ™ Smith ct al. (2008); ™ Peterson et 
al. (2004); Markowitz (2009); Bian & Zhao (2003a); Wang & Zhang (2007); Wandel & Mushotzky (1986); Onken et al. (2003). 



Unifying X-ray winds in Seyfert galaxies 5 



constrained between the inner and outer conical surfaces. 
The assumptions are that the thickness of the wind between 
the two conical surfaces is constant with 5 and that this is 
smaller than the distance to the source. However, as already 
noted in Paper III, we obtain equivalent results including a 
dumpiness factor of ~ Ar/r along the line of sight in the 
spherical approximation case and using a covering fraction 
Cf ~ 0.2/(5, (j>) — 0.4, which is consistent with observations 
(e.g., Blustin et al. 2005; McKernan et al. 2007). Moreover, 
this expression has the same parameter dependencies as that 
recently employed by Crenshaw & Kraemer (2012). Equa- 
tion (3) is actually more conservative, yielding mass outflow 
rates that are roughly a factor of two lower. 

Neglecting additional acceleration of the outflow, i.e. as- 
suming that it has reached a constant terminal velocity, the 
kinetic (or mechanical) power can consequently be derived 
as: 



Ek 



1 



2 



(4) 



The estimates of these parameters are reported in Table 1. 
We also calculated the outflow momentum rate as Pout = 
Moutiiout and subsequently compared it to the momentum 
flux of the radiation field, Prad = Lboi/c. The bolometric lu- 
minosity Lboi is estimated as Lboi = fcboiiion, where fcboi— 10 
(McKernan et al. 2007; Vasudevan & Fabian 2009; Lusso et 
al. 2010). 



4 CORRELATION ANALYSIS 

4.1 Measured absorber parameters 

Initially, we consider the data from all WAs, UFOs and non- 
UFOs together. We compare their ionization parameter, ve- 
locity and column density and search for possible correla- 
tions. The large dynamic range of ~6 orders of magnitude in 
ionization, ~4 in column density and ~3 in outflow velocity 
allows, for the first time, to investigate the general relations 
and trends among these parameters. We fit the power-law 
model log y — a log a; -I- 6 to the datasets using the BCES 
(bivariate correlated errors and intrinsic scatter) regression 
method (Akritas & Bershady 1996) which takes into account 
measurement errors in both the "X" and "Y" coordinates 
as well as the intrinsic scatter. This method has been widely 
used in fitting datasets in the astronomical community (e.g., 
Sani et al. 2011). No a priori dependent variable in the fit- 
ting is assumed and we treat the variables symmetrically. 
Uncertainties on the parameters derived from the fits are 
estimated after carrying out 10,000 bootstrap resamples of 
the data. In order to estimate the significance of these cor- 
relations we calculated the Pearson coefficient Pp. 

The results of the fits are listed in Table 2 and the scat- 
ter plots and best-fit relations are shown in Fig. iQ. All three 
correlations have a high statistical significance(> 6cr). We 
note that the wide dynamic range covered by the observed 
absorption components allows us to fill the whole parameter 
space, with the WAs and the UFOs at the two sides. 

^ In calculating the correlations we used only the points con- 
strained within their errors of measure. Instead, in the figures we 
plot also those with upper or lower limits for completeness. 




-1 



12 3 4 
log^ (erg s^' cm) 



5 6 





5.5 




5.0 




4.5 


'l/l 




E 


4.0: 








3.5 






M 







3.0: 




2.5 




2.0 



(c) 




20.0 20.5 21.0 21.5 22.0 22.5 23.0 23.5 24 

logJVff (cm-2) 

Figure 1. Correlations for measured outflow parameters. Scatter 
plots of log^ vs. logNn (panel a), log^ vs. logtiout (panel b) and 
logNjj vs. logi)out (panel c) for the WAs (red), non-UFOs (green) 
and UFOs (blue). The arrows indicate the lower or upper limits. 
The solid lines represent the best-fit linear correlations and the 
gray shadowed areas indicate the 2a confidence bands. 
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Table 2. Results of the linear regression and partial correlation analysis for the measured absorbers parameters. 



X 


y 


a 


dov(a) 


b 


dev(6) 


scatt 




dof 






z 






^part 




All absorbers together 


log? 


logNn 


0.72 


0.12 


20.00 


0.33 


0.71 


0.73 


55 


1.0 X 10" 


-10 


log-Uout 


0.5 


0.09 


1.0 X 10"" 


10 


log? 


log-Uout 


0.65 


0.10 


6.61 


0.35 


0.52 


0.79 


45 


4.0 X 10" 


-11 


logNn 


0.3 


0.07 


1.0 X 10" 


-4 


logNn 


logfout 


0.69 


0.09 


-6.84 


1.97 


0.45 


0.77 


65 


3.0 X 10- 


-14 


log? 


0.3 


0.09 


5.0 X 10" 


-4 


WAs 


log? 


logNn 


0.73 


0.13 


19.92 


0.29 


0.83 


0.50 




1.5 X 10" 


-3 












log? 


logVout 


0.31 


0.08 


7.19 


0.18 


0.30 


0.69 




1.0 X 10" 


-3 












logNn 


logS^out 


0.48 


0.21 


-2.46 


4.47 


0.30 


0.52 




2.0 X 10' 


-2 






















UFOs and non- 


■UFOs 


















log? 


logNn 


0.62 


0.53 


20.62 


1.98 


0.27 


0.41 




9.0 X 10" 


-2 












log? 


logVout 


0.63 


1.09 


6.72 


4.46 


0.60 


0.23 




2.6 X 10" 


-1 












logNn 


logS^out 


1.43 


1.30 


-23.74 


30.10 


0.79 


0.29 




2.9 X 10" 


-1 













Notes, a and b arc the slope and intercept of the linear correlation fits with standard deviations dev{a) and dcv(fe), respectively, scatt 
represents the internal scatter in units of dex. J?p is the Pearson coefficient, dof is the number of degrees of freedom. -Pnull ^^e 
probability of the null hypothesis that there is no correlation between x and y. z is the third variable against which the partial correlation 
analysis is performed. Tj^ and ap^ are the Kendall's partial correlation coefficient and tlie variance, respectively. Ppart is tire null 
hypothesis probability of the partial correlation. 



In particular, the slope 0.72 ± 0.12 of the correlation 
between the ionization parameter and column density in 
Fig. la suggests that the column density is higher for more 
ionized absorber^ From Fig. lb we observe a linear relation 
with slope 0.65 ± 0.10 between the ionization and velocity, 
indicating that the faster outflows are also more ionized. 
Finally, Fig. Ic shows a relation between column density 
and velocity with slope 0.69 ±0.09, which indicates that the 
faster outflows have also higher columns. In general, these 
correlations suggest that the faster the outflow, the higher 
the column density and ionization parameter. The possible 
systematics and selection effects affecting these relations are 
discussed later in §4.3. 

In order to test whether the correlations between two 
parameters are driven by their mutual dependence on a third 
one, we then performed also a partial correlation analysis 
(Akritas & Siebert 1996). We quantify this effect using the 
Kendall's partial correlation coefficient, tk, which takes into 
account also censored dat4f|. The results of the partial corre- 
lation analysis are reported in the last three columns of Ta- 
ble 2. This test indicates that the correlations between the 
different parameters are only marginally interdependent. 

We performed also a correlation analysis considering 
separately the WAs and the other absorbers (UFOs and 
non-UFOs). As already noted in §2, the distinction between 
UFOs and non-UFOs is only arbitrary and based on their ve- 
locity being higher or lower than 10,000 km/s. In fact, their 
parameter distributions reported in Paper II and Paper III 
are not bimodal. This test was done in order to check if there 

^ In Fig. la we note a possible WA outlier with high column 
density, A'^n ~ 10^** cm"'^, and a relatively high ionization, 
log? ~2.5 erg s~^ cm. As reported in Table 1, this corresponds to 
the Seyfert 1.5 galaxy IRAS 05078+1626. However, given that it 
has not a velocity estimate, this is not considered anymore in all 
the successive relations. 

^ The partial correlation analysis takes into account only the 
points constrained within their errors of measure and those with 
upper limits. 



are significant differences when considering the two popula- 
tions separately. The results are reported in Table 2. Due 
to the large error bars, smaller number of data points and 
reduced dynamic range of the two separated populations, 
their parameter values are less constrained and the signifi- 
cance is lower, especially for the highly ionized absorbers. 
However, a possible difference can be noted in the best- 
fit value of the slope of the relations between log^-logWout 
and logA'n-logWout , the one of the highly ionized absorbers 
being steeper, although they are still consistent given the 
large uncertainties. This point will be addressed later in the 
discussion section. Given the limited quality of the current 
datasets we are not in a position to strongly constrain pos- 
sible differences between the separated correlations. 

Finally, we note that Blustin et al. (2005) attempted a 
similar correlation analysis of the WAs collecting their pa- 
rameters from the literature. They considered a sample of 23 
Seyferts and quasars, finding useful information for only 14 
of them. They found a significant correlation between log^ 
and logAi'n, as reported also by other authors (Holczer et 
al. 2007; Behar 2009). Excluding instead the few fast and 
highly ionized outflows reported at that time, they found 
only a very marginal correlation between log^ and logVout 
and no signiflcant correlation between logUout and logA'n- 
However, it is important to bear in mind that the inability to 
find a significant correlation does not necessarily exclude its 
existence. In fact, after more than seven years of additional 
observations and analysis, we can now deflne a larger sample 
of sources (35) with higher S /N observations, which allowed 
us to increase the number and quality of the reported WAs 
components and, thus, better constrain the correlations be- 
tween the parameters. The inclusion of the UFOs also allows 
to significanly expand the parameter space. 



4.2 Derived outflow parameters 

We checked for possible correlations/trends among the de- 
rived outflow parameters of the WAs listed in Table 1 and 
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Figure 2. Comparison between tlie outflow velocity logUout, the outflow kinetic power log(EK/iEdd) (panel a) and the mass outflow 
rate log(Mout/MEdd) (panel b), respectively. The points relative to the WAs (red open circles), non-UFOs (green filled triangles) and 
UFOs (blue filled circles) are reported. The error bars indicate the upper and lower limits and the points are the average between the 
two. The solid lines indicate the best-fit linear regression curves. 



those of the UFOs and non-UFOs reported in Paper III. 
Given that we can only estimate upper and lower limits, in 
order to carry out the fits we performed 10,000 Monte Carlo 
simulations for each dataset, considering a random value be- 
tween the lower and upper limits for each point (assuming 
a uniform distribution in log space). We then calculated the 
linear regression and the corresponding statistics. 

The average black hole mass of the sources in the 
sample is log{MBH /Mq)~7.4: with a standard deviation of 
~0.7 dex. The difference is not large, but in order to take into 
account the expected scaling of the outflow parameters with 
mass, we normalized the distance to the Schwarzschild ra- 
dius rs = 2GMbh /c^ , the kinetic power to the Eddington lu- 
minosity LEdd = ^nGMsHmpc/aT ~ 1.26x10^* (i\fBH/A/0) 
erg s~^, the momentum to the Eddington momentum rate 
Pgdd = i/Edd/c and the mass outflow rate to the Eddington 
rate MEdd = V^'Edd/c^ , where r; ~ 0.1. 

The parameters and significance of the correlations are 
listed in Table 3. The plots are reported in Fig. 2-4l|- These 
show that, even if the difference between the lower and upper 
limits can be large in some cases, the wide dynamic range of 
the parameters still allows us to estimate significant corre- 
lations/trends among them. The plots of the velocity with 
respect to the kinetic power and mass outflow rate are re- 
ported in Fig. 2. In Fig. 2a we note a strong trend of increas- 
ing the kinetic power for increasing velocity, with a positive 
slope of 1.77±0.14. This is close to the trend expected for 
a roughly constant mass outflow rate. In fact, from Fig. 2b 
we note that, besides the large uncertainties, Mout does not 
vary much with Vout , although there is a possible weak trend 
of decreasing mass outflow rate with increasing velocity with 
slope -0.21±0.14 (see Table 3). 

Instead, in Fig. 3 we show the plots of the different pa- 
rameters with respect to the line of sight projected distance 
in units of Ts . From panels a, b and c we note very significant 
trends of decreasing the ionization parameter, column den- 



sity and velocity with distance, with slopes of — 0.58±0.04, 
-0.44±0.04 and -0.40±0.03, respectively (see Table 3). In 
particular, we checked that the correlation between the ve- 
locity and distance is not an induced relation from the fact 
that the lower limits have been estimated using equation 
(2) in §3, i.e. assuming that the observed velocity is equiva- 
lent to the escape velocity, as this relation is independently 
confirmed using the upper limits alone derived from equa- 
tion (1), as discussed also later in §4.3. In panel d and e we 
can observe a weak increase/decrease of the mass outflow 
rate/momentum rate for increasing distance, with slopes of 
0.16±0.07 and — 0.22±0.07, respectively. Instead, in panel 
f we can note a more pronounced and significant trend of 
increasing the observed outflow mechanical power with de- 
creasing distance, going from the WAs to the UFOs, with a 
slope of -0.60±0.09. 

It should be noted that the whole parameter space 
is essentially uniformly filled, with distances ranging from 
~100 rs from the black hole up to ~kpc scales. The UFOs 
occupy the lower end of the distribution at the smaller dis- 
tances, where the ionization, column and velocity are higher. 
In particular, extrapolating the relations reported in Table 3 
and Fig. 3 to the innermost possible radii of the order of 
log(r/rs) — 1 we obtain that the ionization of the gas reaches 
very high values of log^~5-6 erg s"'^ cm, the column be- 
comes mildly Compton-thick, A^'h — lO^'* cm~^, and the out- 
flow velocity approaches significant fractions of the speed of 
light. The values of the parameters gradually decrease with 
increasing distance, going from the UFOs to the WAs. 

Finally, Fig. 4 shows a significant (99% confidence 
level) linear relation between the radiation momentum rate, 
frad=iboi/c, and the outflow momentum rate of the UFOs, 
Pout, with a slope of 0.76±0.19. The possible physical im- 
plications of this and the previous relations will be briefly 
addressed in discussion section 55. 



* For each of these relations we considered only the points con- 
strained between the errors or upper/lower limits in both the X 
and Y axes. Therefore, the number of points can slightly differ 
from one plot to the other. 



4.3 Possible systematics and selection effects 

In the calculation of the correlations in the previous sections 
we took into account the uncertainties in the ionization pa- 
rameter, column density and outflow velocity as reported 
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Table 3. Linear regression results for the derived outflow parameters. 



X 


y 


a 


dev(a) 


& 


dev(6) 


Rp 


-fnull 




loglJout 


log(EK/I'Edd) 


1.77 


0.14 


-18.89 


1.19 


0.89 


3.8 X 10" 


15 


logtlout 


log(Mout/MBdd) 


-0.21 


0.14 


1.24 


1.19 


-0.22 


1.6 X 10" 


-1 


log(r/rs) 


log? 


-0.58 


0.04 


5.80 


0.21 


-0.85 


8.0 X 10- 


16 


log(r-/rB) 


logATn 


-0.44 


0.04 


24.21 


0.20 


-0.86 


1.3 X 10" 


13 


log(r/rs) 


loglJout 


-0.40 


0.03 


10.47 


0.11 


-0.89 


1.6 X 10" 


16 


log(r/rs) 


log(Mout/MEdd) 


0.16 


0.07 


-1.29 


0.32 


0.38 


5.0 X 10- 


-3 


log(r/rs) 


log(-Pout/PEdd) 


-0.22 


0.07 


-0.44 


0.33 


-0.45 


7.2 X 10" 


-4 


log(r/rs) 


log(-BK/-LEdd) 


-0.60 


0.09 


-0.89 


0.38 


-0.73 


3.1 X 10" 


10 


logPrad 


log-Pout 


0.76 


0.19 


8.68 


6.41 


0.56 


1.0 X 10- 


-2 



Notes, a and 6 are the slope and intercept of the linear correlation fits with standard deviations dev(a) 
and dev(6), respectively, fip is the Pearson coeflicient. fnull probability of the null hypothesis. 



in Table 2. However, when considering such a large datasct 
and especially when collecting results from the literature, it 
is important to bear in mind the existence of possible sys- 
tematics and selection effects. 

As already discussed in Paper II, different assumptions 
of the velocity broadening of the lines can generate some 
variations in the estimated column density. Regarding the 
UFOs, in Paper II we already took into account this effect 
testing different velocity broadciiings for the lines that were 
not resolved and included the larger error bars. For the col- 
umn densities of the WAs collected from the literature, this 
effect is marginal given that the estimates were derived us- 
ing high energy resolution data. Another parameter that can 
affect the estimate of the column density is the assumed el- 
emental abundance. However, as already discussed in Paper 
III, even allowing for a factor of two difference with respect 
to the standard Solar values, the discrepancy in the column 
density is <0.2 dex, within the typical errors of measure. 
The fact that the column densities do not exceed the value 
of A^H— 10^'' cm~^, especially noticeable in Fig. Ic, could in 
principle be affected by the fact that the photo-ionization 
code Xstar can not treat Compton-thick absorbers. How- 
ever, the data do not seem to require significantly higher 
columns as a good spectral modelling of the highly ionized 
UFOs in Paper II was already obtained with columns in the 
range Ah=10^^-10^'' cm"^. 

It is known that the ionization parameter log^ has a 
dependence on the assumed ionizing continuum. In Paper 
III we estimated that the possible uncertainty on the con- 
tinuum slope may induce a maximum systematic error of 
~0.4 dex on the ionization parameter, within the typical in- 
ternal scatter of the relations shown in Fig. 1 (panels a and 
b). Moreover, the large range in observed ionization states 
reduces the importance of this effect when performing the 
linear regression fits. 

Considering the column density and ionization in 
Fig. la, there might be some possible selection effects at 
the two ends of the distribution. In fact, the limited S/N of 
the observations could have hampered the detection of the 
weak spectral features from absorbers with high ionization 
and low column density. On the other hand, the fact that we 
can observe WAs only for type 1 sources could have limited 
the relevance of lowly ionized/neutral absorbers with high 



column densities. However, these latter types of absorbers, 
usually found only in type 2 sources, have intrinsic veloc- 
ity consistent with zero and probably have a different origin 
than the AGN outfiows studied here, such as the --^pc scale 
molecular torus or large scale dust lanes in the host galaxy 
itself. The combination of these two possible selection effects 
could have induced a slight steepening of the log^-logA^H re- 
lation and of the radial density profile successively discussed 
in §5.2. Moreover, these estimates do not take into account 
the possible presence of additional fully ionized material, 
which does not imprint any observable spectral absorption 
features. 

We note that another important parameter of these out- 
fiows is their inclination with respect to the line of sight. 
Unfortunately, the estimate of this parameter is not well 
constrained for each source yet, but the typical inclination 
of type 1 Seyferts is ~30°, with a range between ~10°-60° 
(Wu & Han 2001). Therefore, this ensures that the difference 
among the sources of our sample is not large. However, this 
might contribute to some of the internal scatter observed in 
Fig. 1. Anyway, given the large dynamic range, small differ- 
ences between sources do not significantly affect the derived 
scale relations. 

As already discussed in Paper HI, the expression for the 
mass outflow rate used in §3 introduces a possible system- 
atic source of uncertainty from the assumed typical inclina- 
tion and opening angle of the wind, which however is con- 
strained within a factor of ~0.3 dex (Krongold et al. 2007). 
Instead, the typical systematic uncertainty on the black hole 
masses derived using the reverberation mapping technique 
is <0.5 dex (Peterson et al. 2004). 

Regarding the outflow velocity of the absorbers, as al- 
ready discussed in Paper II, we might be losing some of the 
components with the highest blue-shift due to the fact that 
they are usually detected in the Fe K band at E>7 keV 
where both the energy resolution and effective area of the 
EPlC-pn instrument on board XMM-Newton arc worse and 
the detector has essentially no sensitivity above E~10 keV. 

We also checked that the methods used to estimate the 
the lower and upper limits of the distance of the absorbers 
reported in equations (1) and (2) do not introduce signifi- 
cant systematics in the relations reported in Fig. 2 and 3. 
For instance, the same dependence between the velocity and 
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Figure 3. Comparison between the line of sight projected distance log(r/rs) (in units of the Schwarzschild radius) and the ionization pa- 
rameter log^ (panel a), the column density logA''H (panel b), the outflow velocity logi)out (panel c), the mass outflow rate log(Mout/A'/Edd) 
(panel d), the outflow momentum rate log ( Pout /ffedd) (panel e) and the outflow kinetic power log(ijK/^Edd) (panel f), respectively. 
The points relative to the WAs (red open circles), non-UFOs (green filled triangles) and UFOs (blue filled circles) are reported. The 
solid lines indicate the best-fit linear regression curves. Assuming the typical black hole mass of Mbh — IO^Mq, the distance scale can 
be easily converted from Schwarzschild radii to parsec considering that Ipc ~ 10® rg. 



distance with a slope consistent with ~— 0.4 shown in Fig. 3c 
is independently found also when considering only the upper 
limits of the distance derived from the ionization parameter 
in equation (1). 

All these sources of uncertainties can contribute to the 
significant internal scatter of >0.5 dex observable in the 
plots in Fig. 1. Moreover, we should bare in mind also the 
possibility that inhomogeneities and variability of the ab- 
sorbers and also the distinct analysis methods employed by 
the different authors could contribute to the observed scatter 



as well. A direct, more homogeneous and systematic spectral 
analysis of the WAs, similar to what we have done for the 
UFOs in Paper I and Paper II, would reduce the importance 
of these possible systematics, but this is beyond the scope of 
the present paper. Here we focus more on the global picture, 
noting also that very detailed spectral analysis of WAs could 
be performed only on a handful of sources as they require 
very large and high quality datasets. 

Finally, here we tested for simple linear relations in log- 
space between the absorber parameters, such as ionization. 
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column density, outflow velocity, distance and energetics. 
This represents a good 1st order approximation but there 
might be slight variations in the slope, particularly at the 
low and high ends of the distributions, which could indicate 
different regimes. This possibility will be briefly considered 
later in the discussion section. 



5 DISCUSSION 

5.1 Unification of the X-ray outflows 

In all the tests considered and in Figures 1-3 we find that 
the WAs, non-UFOs and UFOs, within the observational 
uncertainties, show similar relations between their param- 
eters. In particular, the WAs and UFOs lie always at the 
two ends of the distributions, with the intermediate non- 
UFOs in between, and they roughly uniformly cover the 
whole parameter space. When considered all together, wc 
find very significant correlations between their parameters, 
such as ionization, column density, velocity and distance. 
These results strongly suggest that these absorbers, some- 
times considered of different type, could actually represent 
parts of a single general stratified outfiow observed at dif- 
ferent locations along the line of sight. If they were unre- 
lated, the points relative to the different absorbers would 
be mixed and would not display any significant correlations, 
especially when considering several sources together. A sim- 
ple schematic diagram of a stratified accretion disc wind is 
shown in Fig. 5. 

The stratification along the line of sight can easily ex- 
plain the relations in Fig. 1-3. For instance, the fact that 
the faster absorbers are also more ionized suggests that we 
see components that are ejected closer to the black hole. 
The increase in column density for higher ionization can be 
explained with a negative gradient with distance, as also dis- 
cussed in next section when estimating the density profile. 

The observed correlation between the velocity and the 
distance shown in Fig. 3c and reported in Table 3 is con- 
sistent with an approximate slope of —0.40 ± 0.03. This is 
essentially equivalent to the one expected for a wind that 
has an outflow velocity proportional to the escape velocity 
at a given location and, therefore, that is able to escape 
the system and not fall back. Thus, Vout = fvVesc, where 
{2GMIrf^ is the escape velocity and /v ^ 1. There- 
fore, such type of winds are expected to roughly follow a 
relation Vout oc r~^^^ . 

This relation is satisfied in two circumstances: (i) the 
wind observed at each radius was launched relatively close 
to that radius with a velocity proportional the local escape 
speed or (ii) the wind was launched at small radii with a 
velocity close to the escape speed and then decelerated under 
the sole effect of gravity, i.e. ballistically. This latter case is 
probably unphysical because the wind is not expanding in 
pure vacuum, but there is always some interstellar medium. 

In this regard, the observed relation excludes the sce- 
nario of a wind being launched from small radii with a 
constant terminal velocity, which would lead to a flat pro- 
file with radius. However, this last hypothesis is, at least 
partially, ruled out unless some decelerating processes arc 
present, such as shocks or entrainment of some of the sur- 
rounding material. Actually, the fact that the observed pro- 
file is ~— 0.4 instead of exactly —0.5 suggests that some form 



of interaction with the external medium might be present at 
^pc scales. This is supported by the slight increase of the 
mass outfiow rate in Fig. 3d at large distances of >100pc, 
which might indicate some entrainment of surrounding ma- 
terial. This point will be addressed also in §5.4. 

As already discussed in Paper III, the location and char- 
acteristics of the UFOs are indeed consistent with accretion 
disc winds/outflows and the possible direct connection with 
the other non-UFOs and WAs reported here indicates that 
actually all of them could bo identified with the same global 
outflow observed at different locations along the line of sight, 
the WAs being ejected at larger distances, of the order of 
the outer disc and/or torus, the latter possibly representing 
a natural large-scale extension of the disc itself, (see Fig. 5). 
In fact, the correlation analysis and the plots in Figures 1-3 
points toward an underlying connection among these ionized 
absorbers. 

It will be discussed later in §5.3.2 that case (i), in which 
the wind is launched on a wide range of radii with a velocity 
profile proportional to the escape velocity, is probably the 
preferred interpretation and it is directly predicted by strat- 
ified MHD accretion disc wind models. This interpretation 
is valid at least for distances <IClOOpc scales, as additional 
effects of interaction with the surrounding medium might be 
important beyond that (see §5.4). 

Regarding the WAs, we note that the fact that the line 
of sight projected location of some of them is of the or- 
der of the putative obscuring torus at the base of the type 
1/type 2 unification models does not necessarily mean that 
they are produced there, for instance as inferred by Blustin 
et al. (2005). The putative torus is also a large, equatorial 
structure which is difficult to relate with winds that are ob- 
served at relatively small inclinations along the line of sight 
(~30°), as for the Seyfert 1 galaxies discussed here (Wu & 
Han 2001). We also note that some authors even suggested 
that the torus itself could actually be identified with the 
slower and outer regions of a global stratified wind. At these 
locations, dust would be present at the disc surface and it 
could be uplifted and become embedded in the outflow, ex- 
plaining its presence in some WA observations (Konigl & 
Kartje 1994; Elvis 2000; Kazanas et al. 2012). 

Finally, we note that the presence of a photo-ionized 
outflow extending from the inner regions around the black 
hole up to ~kpc scales could be directly related with the 
ionization cones observed through emission lines and images 
in several Seyfert galaxies, some of them being also part of 
this sample (e.g. Storchi-Bergmann et al. 2009; Crenshaw & 
Kraemer 2000; Dadina et al. 2010; Wang et al. 2011a; Paggi 
et al. 2012). 

5.2 Outflow density profile 

An important quantity describing these ionized outflows is 
their radial density profile. Considering all these absorbers 
as different representations of the same phenomenon, we can 
derive an estimate of their average global density profile of 
the form n(r) oc r~". 

From the relation between the column density and dis- 
tance shown in Fig. 3b and reported in Table 3 we have 
ATh = lQ\r/r^Y cm"^ where a = -0.44 ± 0.04 and 
6 = 24.21 ± 0.20. The column density can be expressed also 
as Nyi = n(r)Ar ~ n{r)rs{r/rs), with Ar ~ r. Combining 
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these two equations wc obtain a rough expression for the 
density profile as n{r) ~ no(r-/rs)"~^, with no = lO'/rg ~ 
^^q24.21^3 ^ io5)(Mbh/Mo)-^ 10^°Mg-i cm-^ repre- 

senting the density normalization and Ms = Mbh/IO^Mq. 
Therefore, substituting the observed value of a ~ —0.4, wc 
derive that the density profile has a slope a = 1 — a ~ 1.4. 
Consequently, for a typical black hole mass of ~10'^ Mq and 
for an inner radius of log(r/rs)~l, we obtain a density at 
the base of the outflow of ^10^^° cm~^. 

For comparison, from the relation between the column 
density and the ionization parameter shown in Fig. la and 
following Holczer et al. (2007) and Behar (2009), an esti- 
mate of the radial density profile can be determined us- 
ing also the absorption measure distribution, defined as 
AMD = dNH/d(log^). The AMD is the absorption analog 
of the emission measure distribution (EMD), widely used 
in the analysis of the emission line spectra, and it repre- 
sents the distribution of hydrogen column density along the 
line of sight as a function of ionization. Given the relation 
logNn = alogC -|- & in Table 2 with a = 0.72 ± 0.12, this 
can then be rewritten as AMD = (lO^'alnlO)^" oc Then, 
the slope of the radial density profile can be estimated as 
a — {I + 2a)/(l + a) with uncertainty Aq = Aa/(1 + a)'^ 
(Behar 2009). Substituting the observed quantities in Ta- 
ble 2 we obtain a = 1.42 ± 0.04, consolidating further this 
important result. 

This value is slightly higher than those reported in a 
more detailed analysis of the WAs in a sample of 5 Seyfert 
galajcies by Behar (2009), who nevertheless suggested that 
a slight increase could be present for high ionizations. How- 
ever, it should also be noted that the possible selection ef- 
fects for the absorbers with very low and high ionization and 
column densities previously discussed in §4.3 could induce a 
slight steepening of the density profile estimated in this way. 
Moreover, these estimates do not take into account the possi- 
ble presence of additional fully ionized material, which does 
not imprint any observable spectral absorption features. 

Therefore, in both cases the density profile is n oc r~^'^. 
As already noted by Behar (2009), this scaling rules out two 
simple scenarios, a constant density flow (n oc r°) and, on 
less stringent grounds, a steady, mass conserving spherical 
symmetric radial flow similar to a stellar wind, in which 
the mass outflow rate, the opening angle and the wind ve- 
locity are all constant (n oc r~^). This density profile sug- 
gests that the outflowing gas is more consistent with a con- 
ical/paraboloidal shaped wind instead of a simple spherical 
shell. 

5.3 Acceleration mechanisms 

Once we have established the fact that the ionized absorbers 
can be unified as part of a single, large scale outflow, then 
a fundamental question follows: what is (are) the main ac- 
celeration mechanism(s)? The limited dynamic range in lu- 
minosity for the sources considered here hampers a detailed 
study of possible correlations between the bolometric lu- 
minosity and the parameters of the outflows. In fact, the 
average bolometric luminosity is logLboi— 44.5 erg s~^ with 
a dispersion of only ~0.8 dex. Moreover, the existence of 
different absorber components for each source might render 
the search for such correlations nontrivial. However, we can 
perform some tests. 



In order to get some insights into the possible acceler- 
ation mechanism(s) of the outflow, the velocity alone is not 
a good parameter to compare with the source luminosity, 
while the momentum rate and the kinetic power are better 
because they also take into account the mass outflow rate, 
i.e. the energetics involved. 

Considering the Eddington ratios of the Seyferts in our 
sample showing X-ray outflows we obtain an average value of 
~0.15 with a dispersion of ~0.6 dex. This interval is too nar- 
row to investigate for possible correlations/trends between 
the Eddington ratio and the other parameters. However, we 
note that for a given ratio, the UFOs are always more pow- 
erful than the other outflows. 

More insights might be derived comparing the outflow 
and radiation momentum rates. The outflow momentum 
rate can be expressed as: 

-Pout = Moutv ~ Andmpur'^v'^, (5) 

where Cf is the covering fraction and can assume values be- 
tween and 1. Instead, the momentum rate of the radiation 
field was defined in §3 as Prad ~ I/boi/c. As already intro- 
duced in §4.2 and reported in Fig. 4 and Table 3, there is 
a significant roughly linear (slope a = 0.76 ± 0.19) correlar 
tion between the radiation momentum rate and the outflow 
momentum rate of the UFOs, therefore Pout ~ Pmd- In par- 
ticular, the average value of the ratio between these two pa- 
rameters is consistent with unity, < Pout/P-ad >= 1.6± 1.1, 
indicating a direct proportionality between these two quan- 
tities. 

We checked for the existence of a similar relation also 
for the WAs but we could not constrain it given the large 
errors and scatter. This is due to the fact that their param- 
eters are less homogeneous and they cover a wider interval 
of distances compared to the UFOs. Nevertheless, we can 
estimate a much lower average ratio between their outflow 
and radiation momentum rates of ~0.05, with a standard 
deviation of ~0.26 dex. 

Therefore, the linear correlation observable in Fig. 4 
suggests that there should be a significant exchange of mo- 
mentum between the radiation field and the outfiows or, 
from the relation Mace = L\,oi/ri(P , that the power of the 
outflows is at least related to the accretion rate. 



5.3.1 Radiation pressure 

As we can see from panel a and b of Fig. 3, even if the inner 

part of the flow represented by the UFOs can reach signif- 
icant column densities in the range A'^h — 10^^-10^'* cm^^, 
this material is also highly ionized, with possibly the ma- 
jority of the elements lighter than iron being fully stripped 
of their electrons. Therefore, it is hard, if not impossible, 
to invoke radiation pressure from UV absorption lines as 
the main acceleration mechanism for these flows. This pro- 
cess might instead be more important for winds in bright 
quasars, given their different SED shape and the higher UV 
emission compared to Seyferts (e.g., Proga & Kallman 2004). 

The alternative scenario in which radiation pressure 
could provide significant acceleration to this material even 
if it is highly ionized is then through Compton scattering 
(e.g., King & Pounds 2003; King 2010a, b). In particular. 
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Figure 4. Outflow momentum rate of the UFOs with respect to 
the momentum rate of the radiation field. The error bars indicate 
the upper and lower limits and the points are the average between 
the two. The solid line indicate the best-fit linear regression curve. 



in this case a direct proportionality between the momentum 
rate of the radiation field and that of the outflow would be 
expected: 

Pout.R — CfTcPrad, (6) 

where Tc is the electron optical depth to Compton scatter- 
ing and Cf is the covering fraction of the wind. If this is 
the dominating acceleration mechanism, from the linear re- 
lation between the UFO momentum flux and that of the 
radiation field reported in Fig. 4 we derive that the product 
CfTc should be of the order of unity. 

From the fraction of sources with detected UFOs in Pa- 
per I we derive that the global covering fraction of these 
absorbers is in the range ~0.4-0.6 and therefore we can as- 
sume a typical value of Ct ~ 0.5. Considering the observed 
average column density of the UFOs of A^h — 10^^ cm~^ we 
can estimate an optical depth Tc ~ crTA''H ~0.05. The prod- 
uct of these two values is much lower than ~ unity expected 
from the relation Pout — Prad- However, it should be noted 
that the column density of the UFOs might have been larger 
during the acceleration phase. Extrapolating the relation be- 
tween the column density and the distance shown in Fig. 3b 
and reported in Table 3 we obtain that the column density at 
the innermost radii of log(r/rs)~l is A^h — 10^* cm~^. How- 
ever, this does not take into account the possibile presence 
of some additional fully ionized material that is not visible 
through X-ray spectroscopy. For instance, extrapolating the 
relation in Fig. 3a down to log(r/rs)~l would give rise to 
log^>5 erg s~^ cm. Therefore, we obtain Te ~ 1, indicating 
that the material can be mildly Compton-thick at the base 
of the flow. The product C{Tc is now more close to unity. 

Combining equation (5) and (6) and substituting the 
previously reported definitions for the ionization parameter 
(§2), bolometric luminosity and radiation momentum rate 
(end of §3), we can derive a rough dependence of the outfiow 
velocity with respect to the ionization parameter: 

^ f fcbol V''^ 1/2 f 1/2 /„N 



therefore, for the radiation pressure case we would expect 
'^out oc ri''^^^^^, which is comparable to what found in 
Fig. lb. 

Even if the UFOs have a momentum rate comparable to 
that of the radiation field, it is important to check if the lu- 
minosity of these sources is actually enough to accelerate the 
material to the escape velocity, which is required for such a 
wind to leave the system. Therefore, combining equation (7) 
with the definition of the Eddington luminosity in §4.2 and 
imposing that the wind velocity should be equal or higher 
than the escape velocity at a particular location (see §5.1), 
we obtain: 



where A — Lboi/^Edd is the Eddington ratio, A^h = nAr ~ 
nr is the observed column density and A'^acc is the column 
density of the gas during the acceleration phase. From the 
previously discussed possibility that the column density at 
the base of the fiow (log(r/rs)~l) , where the majority of 
the acceleration should take place, is A'^acc — 10^^ cm^^ 
and considering the average observed column density of the 
UFOs of A'^H — 10^"^ cm~^, we estimate that the radiation is 
capable to accelerate the wind to the escape velocity if the 
Eddington ratio is >0.2. Given that the average Eddington 
ratio of the sources considered here is only A~0.15, we sug- 
gest that radiation pressure might be relevant to accelerate 
the observed outflows. 

The effect of radiation pressure could be increased from 
multiple electron scatterings of the continuum photons if the 
wind opacity would be higher than one or if the luminos- 
ity of these sources would be closer to Eddington. However, 
the maximum opacity at the very base of the wind is only 
Te ~ 1 and these sources are sub-Eddington. Regarding this 
last point, we note that King (2010b) discussed the possi- 
bility that bright AGNs, such as those discussed here, could 
actually be closer to Eddington due to uncertainties on the 
black hole mass and bolometric luminosity estimates. 

Therefore, these evidences indicate that the UFOs may 
be accelerated, at least partially, by the exchange of momen- 
tum with the radiation fleld through Compton scattering. 
This is overall consistent with momentum driven outflow 
models (King & Pounds 2003; King 2010a, b) which actu- 
ally predict the existence of highly ionized outflows with 
velocity ~0.1c and a linear relation between the wind and 
radiation momentum rates. One requirement by these mod- 
els is that the optical depth should be of the order of unity, 
which we find to be possible at the very base of the flow, 
and that the covering fraction should be Cf ~ 1. From ob- 
servations we derive that Cf ~ 0.5 (Paper I), which indicates 
that the outflow has a signiflcant covering fraction and it is 
uncoUimated, nevertheless it is not spherical, but probably 
more consistent with a conical/paraboloidal bipolar wind- 
like shape. 

5.3.2 Magnetohydrodynamic mechanisms 

An additional mechanism that could provide a concur- 
rent acceleration for this highly ionized material is repre- 
sented by magnetohydrodynamic (MHD) processes (Bland- 
ford & Payne 1982; Everett & Ballantyne 2004; Fukumura 
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et al. 2010; Kazanas et al. 2012). We note that MHD mecha- 
nisms are a fundamental requirement for generating the vis- 
cosity in accretion discs (e.g., Abramowicz & Fragile 2011) 
and they are postulated as one of the main heating mecha- 
nisms of the X-ray emitting corona (e.g., Haardt & Maraschi 
1991). Moreover, several Seyfert galaxies do show evidence 
for weak radio jets (Ulvestad & Wilson 1984; Ulvestad et 
al. 2005; Wang et al. 2011b). Therefore, significant mag- 
netic fields might well be present in the inner regions of 
these radio-quiet sources. 

For instance, Fukumura et al. (2010) studied the ion- 
ization structure of self-similar MHD winds off Keplerian 
accretion discs in AGNs. The magnetic field is dragged by 
the rotating disc plasma, and as the wind leaves the disc, 
magnetic torques act on the gas and the wind is magneto- 
centrifugally accelerated. In this case, the wind is found to 
end up with a terminal velocity roughly a few times the ini- 
tial rotational speed. Therefore, a typical feature of these 
winds is that they retain a clear information about their 
launching region. At large distances (>0.1-lpc) the puta- 
tive AGN torus may actually represent an extension of the 
outer accretion disc itself. 

Then, Kazanas et al. (2012) generalized this concept 
and suggested that MHD winds could actually represent the 
fundamental structure at the base of the unified model of 
AGNs. In particular, they derived simple scaling laws for 
these winds and show that they can reproduce several ob- 
served properties of different sources. In particular, the sim- 
ple schematic diagram of a stratified accretion disc wind 
shown in Fig. 5 is very similar to their MHD wind rep- 
resented in their Fig. 8. Magnetocentrifugally accelerated 
winds were also previously studied by other authors, for 
instance Blandford & Payne (1982) and Konigl & Kartje 
(1994), the latter suggesting a similar unification scheme as 
Kazanas et al. (2012). 

When observed with a certain inclination, the line of 
sight intercepts distinct components of this stratified MHD 
accretion disc wind, each with different velocity, column den- 
sity, ionization and projected distance. This model predicts 
that the inner part of the fiow should be faster and more 
ionized, being launched closer to the black hole. The fact 
that the wind is launched from different parts of the disc 
can provide an explanation of the large range in observed 
blue-shifted velocities, ionization and column densities of 
the X-ray absorbers. In this picture, the UFOs and WAs 
would represent the inner and outer part of this accretion 
disc wind, respectively. This is consistent with the relations 
shown in Fig. 1-3 and with the stratified accretion disc wind 
diagram shown in Fig. 5. 

These self-similar solutions provide the simplest means 
of deriving a reasonable MHD wind model and allow to de- 
rive scale relations among different parameters to directly 
compare with observations. A realistic case is definitely more 
complicated than what currently calculated (see Fukumura 
et al. (2010) and Kazanas et al. (2012) for a detailed discus- 
sion of the limitations of these calculations). In particular, 
these models consider the disc only as a boundary condition, 
a fully self-consistent treatment should take into account 
the accretion physics as well. However, they arc sophisti- 
cated enough to provide at least a 1st order characteriza- 
tion of these outflows and allows to investigate the existence 
of underlying relations/trends. Fukumura et al. (2010) fo- 
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Figure 5. Simple schematic diagram of a stratified accretion disc 

wind. The figure is not to scale. This is similar to Fig. 8 of Kazanas 
et al. (2012). The torus, not directly outlined here, may represent 
an extension of the outer accretion disc itself. 



cusod their attention on winds with a radial density profile 
of the type n(r) oc r~™, with a=l. This choice was driven 
by recent observational results on the AMD of some X-ray 
warm absorbers (Behar 2009). Instead, Blandford & Payne 
(1982) adopted a slightly steeper density profile, a = 1.5. 
On the other hand, Konigl & Kartje (1994) also considered 
solutions with radial density slopes of a=l and 1.5, but the 
solution with a=l was preferred because it was representing 
the "minimum energy" configuration. In §5.2 we estimated 
that the general radial density profile of the observed out- 
flows has a slope of roughly a ~ 1.4, which is in between 
these two cases. 

In particular, it is interesting to note the (qualitative) 
resemblance of the scaling relations and dynamic ranges be- 
tween the column density, ionization, velocity and distance 
of the MHD wind model of Fukumura et al. (2010) and 
Kazanas et al. (2012) (shown in the left panel of their Fig. 5 
and Fig. 2, respectively) and our observed relations in Fig. 1 
(panels a and b) and Fig. 3c. For instance, even if they did 
not consider the case of the UFOs in their figure, we note 
the wide range in ionization log^~— 1-5 erg s^^ cm, velocity 
Wout~10-10,000 km s~^ and distance log(r/rs)~4-10. They 
also find the same trends of increasing ionization, column 
density and outflow velocity of tlic> absorbers going from 
large to the small distances to the central black hole. 

Therefore, as already derived for the radiative case in 
§5.3.2, it is interesting to investigate the possible interpreta- 
tion as MHD winds using some more quantitative, although 
still somewhat crude, considerations. The terminal speed of 
magnetocentrifugal winds is proportional to the Keplerian 
velocity at the launching region and it is a few times the 
escape velocity, Wout = (jJVcsc where oj is a factor of the order 
of ~2-3 and it is the ratio between the Alfven radius and 
the wind launching radius (Pudritz et al. 2007). The out- 
flow velocity scales with the Keplerian speed as a function 
of radius, so that the fiow will have an onion-like layering of 
velocities, the fastest inside and the slowest outside. As al- 
ready noted in §5.1, this velocity profile is directly consistent 
with the one observed in Fig. 3c and reported in Table 3. 

It is important to note that, since these winds are ac- 
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celerated by the action of magnetic torques from magnetic 
fields that are embedded in the accretion disc, there is an 
intimate connection between the mass loss rate in the wind 
and the accretion rate onto the black hole. In particular, 
one of the most profound scaling relations for MHD winds 
is represented by the link between the accretion and out- 
fiow rates, Mace ^ w^Mout (Pudritz et al. 2007; Reynolds 
2012). Combining this with the accretion rate expressed as 
Ma.cc ~ Lboi/jyc^, where rj is the radiative efficiency, and the 
definition of the momentum flux for the radiation and the 
wind we can derive the relation: 

^out.MHD — — ^A-ad, (9) 

where /3 = Vout / c, meaning that for a given velocity there is a 
proportionality between the wind and radiation momentum 
rates also for the MHD case. Prom equation (9) and the 
linear relation Pout — P-ad shown in Fig. 4 and discussed 
in §5.3 wc find that the radiative efficiency should be a few 
percent for a typical velocity of the UFOs of /3 ~ 0.1 in 
order to explain this proportionality. Combining equation 
(9) with the expression for the outfiow momentum rate (5) 
and the definition of the ionization parameter we can derive 
a relation between the velocity and the ionization parameter 
for the MHD case as: 

Wout,MHD =^ 2 ( ^2°A ) ^' (10) 

therefore, in this case we would roughly expect a direct pro- 
portionality between the velocity and ionization Wout oc S./v- 
As already done for the radiation pressure case, it is impor- 
tant to check the general conditions under which an MHD 
wind can actually form, imposing that the outflow velocity 
from equation (10) is equal or higher than the escape veloc- 
ity (§5.1). Using the approximation ^ ~ Lion/Nur, the def- 
inition of the Eddington luminosity and the Schwarzschild 
radius, we obtain: 

A > 2V^Te{r)r]0j^Ci. (11) 

where A = Lboi/iEdd is the Eddington ratio, Te{r) is the 

Compton scattering optical depth at the radius f = r/r,. 
Solving this equation for r; and considering an average 
A ~ 0.15, an optical depth at the radius logf~l of ro(f)~0.1- 
1 (see §5.3.1 and Table 3) and the typical values of a;~2-3 
and Cf ~ 0.5, we obtain an estimate of the radiative effi- 
ciency of the accretion disc of the Scyfcrt galaxies consid- 
ered here in order to be able to efficiently generate MHD 
winds: ?? < X/2V^Te{f)oj^C( <0.05-0.1. This value of the 
radiative efficiency is comparable to the typical one derived 
for quasars of 77--O.I (Soltan 1982; Elvis et al. 2002; Davis & 
Laor 2011). Regarding Scyfcrt galaxies in particular, some 
authors suggested that the radiative efficiency for some of 
them could be slightly lower than that (Bian & Zhao 2003b; 
Panessa et al. 2006; Singh et al. 2011). 

Therefore, from these considerations on the wind ener- 
getics and the consistency with the expected overall struc- 
ture/geometry and velocity profile we derive that the ob- 
served outflows could be effectively accelerated by MHD 
processes. 



However, it is important to note that this is only a par- 
tial conclusion because a complete characterization of these 
outflows would require the combined treatment of both ra- 
diation and MHD effects, both important in AGNs. Some 
attempts in this direction have been reported in the litera- 
ture (Proga 2000, 2003; Everett 2005; Ohsuga et al. 2009; 
Reynolds 2012). If these processes are acting simultaneously, 
we could naively expect to observe changes between these 
two regimes. For instance, from equations (7) and (10) we 
can speculate that the velocity should roughly show a pro- 
portionality to 5^^^ from the radiation pressure term and to 
^ from the MHD part. It is tempting to compare this with 
Fig. lb, where we can sec that the residuals of the data 
with respect to the linear fit are possibly consistent with 
a similar change of slope between these two regimes going 
from lower to higher ionization. In fact, the linear regression 
fit is consistent with an intermediate slope of ~0.65. Very 
close to the black hole and for higher ionization/velocities, 
the MHD regime should always be the dominant part and 
it might actually enter a regime that is eventually responsi- 
ble for the acceleration of jets. Further away from the black 
hole and for lower inization/velocities, radiation forces may 
contribute more to the properties and dynamics of the fiow, 
depending of the actual state of the material and also the 
local disc luminosity. 

5.4 Energetics and feedback impact 

Having established that the correlations among the differ- 
ent outflow parameters suggest an interpretation as a strati- 
fied, large-scale wind with most probably both radiation and 
MHD mcchainsms having a role in the acceleration, then the 
next step is to, at least qualitatively, investigate the contri- 
bution of these outflows on the expected AGN feedback. 

An important parameter to quantify the effective feed- 
back contribution of winds/outflows in bright AGNs is the 
ratio between their mechanical power and the bolometric 
luminosity. Extensive numerical simulations demonstrated 
that this value should be around ~5% and can be as low as 
~0.5% (e.g., Di Matteo et al. 2005; Hopkins & Elvis 2010). 
Gaspari et al. (2011a, b, 2012b) recently demonstrated that 
outflows with these powers are indeed adequate to produce 
sufficient feedback to quench cooling flows in elliptical galax- 
ies and to significantly eject gas without overheating the 
gala:xy itself. Since Seyferts are generally hosted in spiral 
galajcies, which typically have less dense environments com- 
pared to ellipticals, then the threshold could be even lower 
than ~0.1% (Gaspari et al. 2011a, b, 2012b). 

Simulations of AGN outfiows with characteristics equiv- 
alent to UFOs have also been independently demonstrated 
to be able to significantly interact not only with the inter- 
stellar medium of the host galaxy but possibly also with the 
intergalactic medium. They can provide a significant contri- 
bution to the quenching of cooling flows and the inflation of 
bubbles/cavities in the intergalactic medium in both galaxy 
clusters (e.g., Sternberg et al. 2007; Gaspari et al. 2011a) 
and esiKxially groups (e.g., Gaspari et al. 2011b). 

In Fig. 6 we plot the kinetic power of the outflows with 
respect to the bolometric luminosity. Wc note again that 
the distribution of points seems rather continuous from the 
WAs to the UFOs, the latter having a higher power for a 
given source luminosity. As discussed in Paper HI and as 
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Figure 6. Outflow kinetic power witii respect to tiie bolometric 
luminosity. Tlie points correspond to the WAs (red open circles), 
non-UFOs (green filled triangles) and UFOs (blue filled circles), 
respectively. The error bars indicate the upper and lower limits 
and the points are the average between the two. The transverse 
lines indicate the ratios between the outfiow mechanical power 
and bolometric luminosity of 100% (solid), 5% (dashed) and 0.5% 
(dotted) , respectively. 

evident from Fig. 6, the UFOs have indeed a mechanical 
power clearly higher than 0.5% of the bolometric luminos- 
ity, with the majority of the values around ~5%. However, 
we should note that, as recently reported by Crenshaw & 
Kraemer (2012), some of the WAs might actually reach the 
~0.5% level as well when their components are co-added 
together. Therefore, these outflows, and in particular the 
UFOs, are clearly the most promising candidates to signif- 
icantly contribute to the AGN feedback besides radio jets 
(e.g., Fabian 2012). 

Theoretically, feedback from AGN outflows has been 
demonstrated to clearly influence the bulge star formation 
and SMBH growth and possibly also to contribute to the 
establishment of the observed SMBH-host galaxy relations, 
such as the Mbh-o" (e.g.. King 2010a; Ostriker et al. 2010; 
Power et al. 2011; Zubovas & King 2012; Faucher-Giguere 
& Quataert 2012). Similar and possibly even more massive 
and/or energetic outflows might have influenced also the for- 
mation of structures and galaxy evolution through feedback 
at higher redshifts, close to the peak of the quasar activity at 
z ~ 2 (Silk & Rees 1998; Scannapieco & Oh 2004; Hopkins 
et al. 2006). 

From Fig. 3a we can note that on ~kpc scales the out- 
flow is very lowly ionized (log^ < erg s~^ cm) and it 
could represent the possible conjunction point with large- 
scale neutral/molecular outflows recently found in some 
sources in other wavebands (e.g. Nesvadba et al. 2008; Sturm 
et al. 2011). Several models have been suggested in order 
to explain their origin, but essentially all of them rely on 
a two step process in which an initial ~sub-pc scale fast 
{vout > 1,000 km s~^) AGN accretion disc outflow per- 
turbs/shocks the interstellar medium, sweeping it up on its 
way and then decelerates/cools (e.g. King 2010a; Zubovas & 
King 2012; Faucher-Giguere & Quataert 2012). Therefore, it 
is tempting to check for possible evidences of these efTects in 
our correlation plots. Their intensity should be more promi- 
nent at large distances from the black hole. Besides the large 
uncertainties, in Fig. 3d we can note a slight increase of the 



mass outflow rate at >100pc (> 10*rs) scales, possibly sug- 
gesting some entrainment of surrounding material by the 
wind. In Fig. 3e and f we can also see a slight increase of 
the wind momentum rate and mechanical power at those 
locations. These evidences are roughly consistent with the 
relations reported in Fig. 4 of Faucher-Giguere & Quataert 
(2012), who performed a detailed study of the interaction 
of AGN winds with the surrounding environment and the 
difTerent regimes of momentum/energy conservation as the 
resulting shocked material propagates to large distances. 

Observationally, we note that evidences for AGN feed- 
back activity driven by outflows /jets improved signiflcantly 
in recent years but there are still significant uncertain- 
ties, especially regarding the link between the observed 
phenomenologies at small (~pc) and large (~kpc) scales. 
Promising results on this line have been recently reported 
for a few Seyfert galaxies, with the detection of bubbles, 
shocks and jet/cloud interaction, some being also part of our 
sample (e.g., Wang et al. 2010; Pounds & Vaughan 2011). 



6 CONCLUSIONS 

In order to investigate the possible relations between the 
ultra- fast outflows (UFOs), mainly detected in the Fe K 
band through Fe XXV/XXVI absorption lines, and the soft 
X-ray warm absorbers (WAs) we performed a literature 
search for papers reporting the analysis of the WAs in the 35 
type 1 Seyferts of the sample defined in Paper I. The main 
results of our study are: 

• The fraction of sources with reported WAs is >60%, 
consistent with previous studies. The fraction of sources with 
UFOs is >34%, > 67% of which showing also WAs. 

• We reported the main observed WA parameters, such 
as ionization, column density and outflow velocity. Then, 
from these values, we estimated also the mass outflow rate, 
momentum rate and mechanical power. 

• The large dynamic range obtained when considering all 
the parameters of these absorbers together allows us, for the 
first time, to estimate significant correlations among them. 
We find that the closer the absorber to the black hole, the 
higher the ionization, column density, velocity and therefore 
the mechanical power. In particular, in the innermost part 
of the flow, at distances of log(r/r-s)~l, we flnd that the 
material can be mildly Compton-thick, A^h ~ 10^^ cm~^, 
highly ionized, log^ ~ 5 erg s~^ cm, and the velocity can 
approach significant fractions of the speed of light. 

• In all the tests the absorber parameters uniformly cover 
the whole parameter space, with the WAs and UFOs lying 
always at the two ends of the distribution. This strongly 
indicates that these absorbers, sometimes considered of dif- 
ferent type, could actually be unified in a single, large scale 
stratified outflow observed at different locations along the 
line of sight. The UFOs are likely launched from the inner 
accretion disc and the WAs at larger distances, such as the 
outer disc and/or torus. See Fig. 5 for a simple schematic 
diagram of such a stratified wind. 

• Given the high ionization and velocity of the outflows, 
and a linear relation between the outflow and radiation mo- 
mentum rates, we argue that the only two viable acceler- 
ation mechanisms are radiation pressure through Compton 
scattering and MHD processes, the latter playing a major 
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role. In particular, the overall structure/geometry is more 
consistent with a stratified MHD wind scenario. 

• Finally, as already discussed in the previous Paper III, 
here we confirm that these outflows, with the UFOs rep- 
resenting the most energetic part, have a sufficiently high 
mechanical power (>0.5% of the bolometric luminosity) to 
provide a significant contribution to AGN feedback. In this 
regard, we find possible evidences for the interaction of the 
AGN wind with the surrounding environment on large-scales 
from the correlation plots. 

In the future, in order to better quantify the predomi- 
nance of radiation pressure or MHD driving, we should ex- 
tend the sample to sources with lower and higher Eddington 
ratios. One would naively expect outflows in lower lumi- 
nosity AGNs to be dominated by MHD processes while for 
bright quasars at the other end of the distribution probably 
radiation pressure is more important. Seyfert galaxies might 
represent an intermediate case between these two, also sug- 
gested by their Eddington ratio of ~0.1. In this respect, we 
plan also to directly test these hypotheses fitting the data 
with detailed radiation and MHD wind models (e.g., Sim et 
al. 2010; Fukumura et al. 2010). 

From an evolutionary point of view, it will be interesting 
to compare the characteristics of similar outflows found in 
higher redshift quasars in X-rays (Chartas et al. 2002, 2003, 
2009; Giustini et al. 2011; Lanzuisi et al. 2012) and also 
powerful, large-scale outflows detected in other wavebands 
(e.g., Nesvadba et al. 2008; Sturm et al. 2011). 

Finally, we anticipate that the unprecedented high en- 
ergy resolution and sensitivity in the wide E~0.1-10 keV 
energy band of the microcalorimeter on board the upcom- 
ing Astro-H mission will provide important improvements 
in this field, allowing for the first time to simultaneously 
study in detail the absorbers in a wide range of ionization 
states, column densities and velocities. We also note that 
the unprecedented effective area of about 10 m^ at 8 keV 
with even moderate energy resolution of ~250 eV, such as 
the one proposed by the LOFT mission to the ESA Cosmic 
Vision, would allow to detect UFOs in bright local AGNs 
with high significance and at velocities up to ~0.7c thanks 
to the extended energy bandpass at higher energies. 
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